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Introduction 


The  aim  of  this  phase  of  the  project  was  to  create  a  tool  for  linking 
the  flow-pictures  of  the  disturbed  two-dimensional  mixing-layer  and  the 
phase- locked  averaged  data  calculated  from  experimentally  measured  velocity 
signals.  We  could  thus  relate  the  quantitative  analysis  of  our  experimen¬ 
tal  data  to  the  visually  observed  phenomena.  Of  particular  interest  is  the 
presence  of  a  pairing  process  of  vortices  in  the  experimental  data. 

There  are  two  theoretical  approaches  to  deal  with  the  mixing-layer 
type  of  flow.  One  is  based  upon  stabilty  theory,  and  the  other  tries  to 
solve  the  basic  flow  equations  by  various  mathematical  models. 

Gorcos  and  Sherman  (1976)  created  a  numerical  model,  and  used  it  to 
examine  the  roll  up  process  of  the  shear  layer  into  a  vortex.  They  pro¬ 
duced  a  series  of  graphs  showing  iso-dynes  (iso-vorticity  contours)  and 
streamlined,  forming  a  "Cat  Eye"  structure,  with  a  vorticity  concentration 
insidd,  connected  to  its  neighbours  by  vorticity  concentrated  braids. 
Acton  (1976)  suggested  a  model,  consisting  of  vortex  elements  which  pro¬ 
duced  the  definition  required  to  observe  the  evolution  and  coalescence  of 
the  large  eddies.  Ashurst  (1976)  employed  the  method  of  discrete  vortices 
in  order  to  simulate  the  two-dimensional  mixing-layer.  He  produced  a  ser¬ 
ies  of  streakline  patterns,  showing  the  time  development  of  the 
aixing-layer . 
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Riley  and  Metcalfe  (1980)  performed  direct  numerical  simulations  of 
the  non-  linear  Navier-Stokes  equations  with  pseudo spectral  numerical  meth¬ 
ods.  They  employed  no  explicit  turbulence  modelling,  and  succeeded  in  mo¬ 
delling  the  vortex-  pairing  mechanism.  They  found  that  the  presence  of  the 
out-of-phase  subhannonic  of  the  most  unstable  linear  eigenmode,  is  essen¬ 
tial  for  the  vortex  pairing  roll-  up  process.  By  eliminating  the  subhar- 
monid,  mixing-layer  growth  stops  and  negative  turbulent  energy  production 
occurs,  along  with  a  countergradient  momentum  flux.  They  suggested  that  in 
the  laboratory  experiments,  the  out-of-  phase  subharmonic  is  suppressed 
significantly,  delaying  vortex  pairing. 

Inoue  (1983)  simulated  turbulent  mixing  layers  by  a  discrete  vortex 
method  obtaining  streaklines  of  fluid  particles  as  well  as  those  of  vor¬ 
tices.  He  found  two  stages  in  the  dev  elopement  of  the  mixing  layer:  a 
stage  of  clustering  where  the  vortices  interact  with  each  other  and  tend  to 
form  separate  clusters  of  discrete  vorticed,  and  a  vortex  pairing  stage. 
He  compared  his  results  with  experiments  and  found  them  to  be  in  good 
agreement  both  qualitatively  and  quantitatively. 

The  other  physical  quantity  to  be  utilized  is  the  stability.  Michalke 
(1964,  1965)  used  this  approach  to  predict  the  temporal  and  spatial  growth 

of  disturbances  in  an  inviscid  two-dimensional  mixing-layer,  approximating 
the  mean- velocity  nature  by  a  hyperbolic- tangent  profile.  His  results  show 
the  roll-up  of  a  vorte*,  both  by  iso-vorticity  contours  and  streaklines. 
Monkevltz  and  Huerre  (1982)  studied  the  influence  of  the  velocity  ratio  on 
the  spatial  Instability  of  mixing-layers,  using  both  hyperbolic-tangent  and 


Blausius  mean-velocity  profiles.  They  found  that  the  maximun  growth-rate 


is  approximately  proportional  to  the  velocity  ratio.  Crighton  and  Gaster 
(1976)  used  the  multiple-scales  expansion  method  to  analyze  the  develope- 
ment  of  a  slowly-diverging  jet  flow.  Gaster,  Kit  and  Wygnanski  (1984)  ap¬ 
plied  the  same  method  to  the  two-dimensional  mixing-layer,  approximating 
the  mean- velocity  behaviour  by  modified  hyperbol ic- tangent  profiles,  and 
showed  an  agreement  between  the  mixing  layer  growth-rate  results  thus  obta¬ 
ined,  and  experimental  results  in  the  weakly-disturbed  turbulent 
mixing-layer.  They  also  showed  the  similarity  between  the  eigenfunction 
profiles,  of  both  the  amplitude  and  phase,  to  the  stream-wise  velocity  com¬ 
ponent. 


In  'this  report  we  will  corapare-^oth 


the  longitudinal  and  lateral  velo¬ 


city  components,  as  well  as  the  vorticity  distribution,  using  the 
multiple-scales  expansion  method.  Simultanous  measurements  of  two  velocity 
components  formed  the  experimental  basis  for  this  comparison. 
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I  Data  Acquisition  and  Numerical  Calculations 

J 


The  first  step  in  the  nunerical  computations  was  to  translate  the  vol- 
tageA,  recorded  on  the  magnetic  tape  during  data  acquisition,  to  the  two 
velocity  components,  using  the  calibration  coefficients.  This  procedure 
produced  the  instantenous  values  of  Q  and  oi  ,  and  the  velocity  components 
U  and  V  were  calculated  from  the  following  relations: 

U=Q  cost  ©(  )  and  V=Q  sin(o<  ) . 

2.1  Phase-locked  averaged  velocity  components 

The  phase-locked  averaged  velocity  vectors  were  calculated  from  the 
sampled  velocity  vectors,  in  order  to  analyze  the  periodic  flow  phenomena, 
having  the  disturbance  frequency,  or  its  first  subharmonic. 

The  nunerical  representation  of  the  phase-locked  averaged  signal  is: 

»  , 

*<U(j)>=  u(®n+(j-1  )At+i«NPP) 

<«/  1  u 

Where  u  is  the  original  vector  containing  the  velocity  values  sampled,  A  t 
is  the  time  Interval  between  data  pointA,  NPP  is  the  number  of  data  points 
in  each  time-period  (T)  of  the  frequency  analyzed  (i.e.  NPP=  T  /  At),  <p  0 
is  an  arbitrary  initial  phase,  and  N  is  the  nunber  of  time-periods  of  the 
analyzed  frequency,  contained  within  the  original  velocity  vector. 


The  determination  of  (f>  ^  and  NPP  required  the  simultaneous  sampling 
and  recording  of  both  the  sine-wave  signal  activating  the  oscillating  flap, 
and  its  first  subharmonic  signal,  produced  by  the  electronic 
frequency-divider  circuit,  designed  to  produce  a  synchronized  square-wave 
signal . 


The  following  sketch  illustrates  the  phased-locked  averaging  method. 


The  initial  phase  points  $  q  were  defined  as  the  cross  points  of  the 
subharmonic  square-wave  signal  ,  having  a  negative  slope,  with  the 
zero-voltage  level,  and  NPP  was  calculated  as  the  mean  nunber  of  data  po¬ 
ints  between  two  adjacent  initial-phase  points. 


The  result  of  an  averaging  process  of  a  signal,  dominated  by  a  certain 
frequency  (a  fundanental  frequency)',  phase-locked  averaging  it  with  respect 
to  its  first  subharmonic,  is  a  signal  having  the  subharmonic's  time-period, 


.jvt.--.  .  r. ^  -  . . ... 


-  **  uJk  '  •£  '  A-  WL  aV  A 
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containing  two  consecutive  signals  of  the  fundamental's  time-period.  All 
high  frequency  interferences  (turbulent  behaviour)  are  filtered  out.  If 
the  subhannonic' s  behaviour  is  phase-locked  to  the  fundamental  frequency, 
the  two  signals  would  be  different  from  one  another  (shown  schematically  as 
[d],  in  the  sketch  above).  A  non- phase-locked  subharmonic  behaviour*,  or 
it's  non-existence  ,  would  produce  identical  fundanental  signals  [  e] . 


The  length  of  each  sampled  velocity  vector  was  512  points  (equivalent 
to  a  period  of  0.25  sec)',  and  100  such  vectors  were  acquired.  Thus ,  the 
subhannonic  phase-locked  average  was  calculated  from  *100  events. 


Other  flow  parameters,  such  as  the  mean  velocity  profiles,  the  R.M.S. 
values  and  the  Reynolds  stresses,  were  calculated  from  the  original  time 
dependent  velocity  vectors,  and  spectral  analysis  was  performed  as  well, 
using  a  FFT  procedure. 


2.2  Reynolds  stress 


Three  different  ways  of  calculating  the  Reynolds  stresses  (or  rather 
the  u-v  product)  were  employed,  each  describing  a  different  aspect  of  these 
stresses: 


a)  The  Reynolds  stress  by  its  definition,  from  the  original  time  de¬ 


pendent  velocity  vectors: 


'£7(X,Y)=-  Np(u(X,Y,t)*v(X,Y,t)) 


WWW 


-r.  _ 


where  u(X, Y,t)=U (X, T,t)-U (X,Y)  and  v(Xf Y,t)=V (X, Y,t)-V ( X, Y ) 


b)  The  phase-locked  averaged  u-v  product,  producing  a  value  containing 

the  influence  of  all  frequencies  present  at  the  phase-locked  averaged  velo¬ 
city  components'  vectors  (frequencies  higher  than  the  phase-locking  fre¬ 
quency  that  were  not  filtered  out  during  the  phase-locked  averaging  pro¬ 
cess):  <  T(x*,  y/)>  -  -  S  <-' v  Cx.,*^o  >•<  v  o<icr 

Where  T  is  the  time-period  of  the  phase-locking  frequency. 

c)  The  Reynolds  stress  spec  trim. 

In  order  to  obtain  both  the  power  spectra  and  the  Reynolds  stress 
spectra,  we  have  used  a  fast  Fourier  transform  operator  (FFT).  The  result 
was  a  vector  containing  both  the  real  and  imaginary  coefficients  of  various 
frequencies,  for  each  velocity  component. 

o£*  (U )  =  C  ^  ( w)  where  C.(j)sA 

£  (V)  =  C2(w)  where  j )  =A2 (  j)+iB2(  j) 

where  A(j)  and  B(j)  are  the  sine's  and  the  cosine's  coefficients  of  the  j*'*1 
frequency,  respectively.  The  sampling  frequency  (2.0*48  KHz),  and  the 

length  of  each  channel’s  sampling  vector  (512  points)',  define  both  the 
lowest  frequency,  and  the  spectra's  resolution  to  be  4  Hz,  and  the  highest 
frequency  to  be  1.02*4  KHz. 


The  power  spectra  was  calculated  by  multiplying  the  coefficients'  vec- 
tor  C^(w)  by  its  complex  conjugate  C^w)  ; 


’l(w):Ci(w)»C» 


(w)  where 


i1(j)  =  A1(j)2t-B1(j)‘ 


a 
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The  Reynolds  stress  spectra  was  achieved  by  calculating  the  real  part 
of  the  product  of  one  velocity  components'  vector  by  the  other's  complex 
conjugate  : 

T(w)  =  0.5  REAL[C1(w)*C*(w)] 
where  T(j)  =  -0.5 [A1 ( j)*A2( j)+B1 ( j)#B2( j) 

2. 3  Vorticity  field 

The  phase-locked  averaged  vorticity  field  was  computec  r  numerical 
differentiation  of  the  phase-locked  averaged  velocity  vectors. 

The  resulting  phase- locked  averaged  vorticity  signal  <w(X,  Y,  T)>was 
used,  both  for  plotting  the  iso-dynes  (iso-vorticity  contours),  and  for 
further  mathematical  operations  (i.e.  Fourier  transformation  in  the  dis¬ 
turbance  frequency  in  order  to  obtain  the  vorticity  eigenfunctions  and  com¬ 
pare  them  with  the  results  of  the  stability  theory). 


2.^  Streaklines 


A  streakline  connects  , by  definition ,  all  fluid  particles  that  have 
passed  through  a  fixed  point  in  space  (a  source  point),  during  a  given  per¬ 
iod  of  time. 


In  order  to  establish  and  draw  a  streakline  from  the  velocity  vectors, 
given  in  the  Eulerian  representation ,  one  must  transfer  it  to  the  Lagrangi- 
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an  representation  ,  i  .e .  calculate  the  particles'  pathlines  given  by: 

r  =  r(R,t) 

where  R  is  the  position  vector  of  the  particle  at  t=0. 

Providing  the  Jacobyan  of  the  equations'  matrix  has  a  non-zero  value 
the  following  relation  can  be  written: 

R=R  (r  ,t) 

A  streakline  at  a  given  time  t  can  be  directly  calculated  by  the  fol¬ 
lowing  expression: 

r  =  r[R(r0,t0)ft'] 

where  rQ  is  the  position  vector  of  a  source  point,  t  is  a  parameter  defin¬ 
ing  the  time  at  which  the  streakline  is  to  be  plotted  and  tQ  is  the  time 
variable  0  <  tQ  <  t  . 


I 

I 

J 

'i 
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Another  way  of  representing  a  pathline  in  a  two-dimensional  flow  field  , 
is : 

U=dX/dt  ;  V  =dY/dt  with  X(t=t0)=X0  and  Y(t=t  )=YQ 

where  (X^.Y^)  are  the  space  coordinates  of  a  source-point. 

We  restricted  our  discussion  to  a  row  of  source-points,  displaced 
vertically  at  the  first  X-station  (X=200  mm).  Thus  we  can  represent  the 
pathlines  as  : 

rl  =  r1(R,tQ,t) 

where  R  is  the  source-points'  position  vector  parameter,  t^  is  the  intial 
time  parameter  and  t  is  the  time  variable.  As  the  flow  field  considered  is 
periodic,  with  a  time-period  T,  it  follows  that: 


- 1  - ’.v 


The  pathlines  were  calculated  by  miner  ical  integration  of  the 
phased-locked  averaged  velocity  vectors. 

Each  source  point  (Xq,Yq)  produced  NPP  pathl ines-the  number  of  data 
points  in  the  phase-locked  averaged  velocity  vector.  The  pathline  vectors 
were  identified  by  an  index  j  defined  as  ^=1+^/  t  where  tQ  is  the  time 
the  particle  passed  through  the  source  point  and  t  is  the  time  interval 
between  data  points. 

Since  At  is  determined  by  the  sampling  frequency  (  t=1/f),  the  in¬ 

tegration  led  us  to  space  coordinates  other  than  our  measuring  stations. 
Interpolation  techniques  were  employed  in  order  to  establish  the  velocity 
values  at  those  points.  A  spline  interpolation  peocedure  was  first  used  in 
the  transverse  direction  (Y-  axis)',  and  a  linear  interpolation  method  was 
then  applied  in  the  stream-wise  direction  (X-axis).  Since  the  time-period 
corresponds  to  NPP  data  points,  all  time  considerations  reduce  to  a  proper 
indexing  of  the  points  in  the  phase-locked  averaged  velocity  vectors. 

The  function  r^  describes  the  position  of  all  particles  that  had 
passed  through  a  source  point.  Writing  : 

r1=rl(R,t0,t*) 

where  R  and  t*  are  now  the  parameters  of  the  source  point  position  and  the 
plotting  time  respectively,  and  t^  is  the  variable,  (  0  <  tQ  <  t  ) ,  we  re¬ 
present  the  desired  streaklines.  The  periodic  flow  field  demands  that: 
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r^R.tg.t*)  =  r1  (R,tp,t#+nT) 


n=  1,  2. 


Thus  the  streaklines'  computation  reduces  to  a  proper  picking-up  procedure 
of  data  points  from  the  pathlines'  vectors. 


If,  for  exanpld,  the  streakline  at  the  time  corresponding  to  the  i 

* 

point  of  the  phase-locked  averaged  vector  is  to  be  plotted  (i.e.  t  =i*  t)  , 
the  following  points  are  to  be  picked  up  from  the  pathline  vector  indexed 


ith,  NPP+ith, 2*NPP+  ith . n*NPP+ith, 

from  the  pathline  indexed  j=2  the  following  points  are  to  be  picked-up: 

(i-1  )th,(NPP+i-1  )th  ,(2*NPP+i-1  )th,  .  .  .  ,  (  n*NPP+i-1  )th  , 
from  the  pathline  indexed  j=j  the  following  points  are  to  be  picked-up: 

(i->1  )th,(NPP+i->1  )th,(2«NPP+i-j+1  )th,  .  .  .  ,(n*NPP+i-j+1  )th  , 
and  from  the  pathline  indexed  j=i: 


1st , K'PPth,  2«NPPth . n«NPPth, 


The  interval  of  NPP  points  between  each  two  points  picked  up  from  the 
sane  pathline  corresponds  to  the  time  period  (T)  of  the  flow. 
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The  Disturbed  Mixing-Layer 


3. 1  Streaklines 

One  of  the  main  objectives  of  this  work  was  to  compare  the  streaklines 
computed  from  experimentally  measured  velocity  signals  with  conventional 
flow  visualization  pictures,  thus  providing  a  link  between  the  qualitative 
observations  of  flow  (flow  visualization)  and  quantitative  measurements 
taking  into  account  the  dynamics  of  the  flow. 

A  collection  of  streaklines,  each  one  originating  from  a  different 
source-point,  all  displaced  vertically  at  the  first  measuring-station 
(X=200  mm),  are  shown  in  Fig.  1.  A  comparison  between  this  streaklines' 
pattern  and  a  smoke  picture  is  presented  in  Fig.  2,  and  the  resemblence 
between  the  two  is  clearly  seen.  The  snoke  pictures  were  taken  using  stro¬ 
boscopic  illumination  at  the  disturbance  frequency,  thus  it  represents  a 
"phase- locked"  averaged  picture  of  36  events. 

The  three  regions  of  growth  described  in  our  first  report,  and  shown 
in  Fig.  3t  can  now  be  discussed  in  terms  of  the  streaklines.  Region  I, 
which  is  marked  by  a  rapid  growth  of  the  shear  layer1,  can  be  characterized 
by  amplification  of  disturbances  in  the  lateral  direction.  Its  length  is 
therefore  equal  to  the  distance  required  for  the  vortex  sheet  to  start  rol¬ 
ling  up  into  lumps.  We  define  a  lump  as  a  structure  in  which  the  streak- 


lines  fold  over  backwards.  The  lateral  dimension  of  the  vortex  lumps  in 
region  II  remains  unchanged  while  their  inclination  to  the  flow  axis  is 
being  altered.  The  thin  vortex  sheet  which  seperates  adjacent  limps,  and 
sometimes  is  referred  to  as  the  "braid",  is  stretched  in  the  process.  In 
region  III  the  vortex  lunps  are  leaning  downstream  on  the  high  velocity 
side  of  the  shear  layer;  their  inclination  angle  to  the  axis  of  the  flow 
continuously  decreases  with  increasing  distance  from  the  splitter  plate. 
The  elongation  of  the  vortex  limps  and  the  change  in  their  inclination 
angle  causes  one  vortex  lunp  to  overlap  another.  The  increase  of  the  aver¬ 
age  momentum  thickness  in  region  I  is  concurrent  with  an  increase  in  the 
lateral  displacement  amplitude  of  a  streakline  originating  at  the  heart  of 
the  mix ing- layer* ,  which  can  be  regarded  as  an  interface  separating  the  two 
mean  streams.  In  region  III  the  duration  of  passage  of  each  vortical 
structure  at  a  given  station  increases  without  a  concommitant  increase  of 
its  lateral  dimension.  This  increase  is  responsible  for  the  widening  of 
the  momentum  thickness.  In  region  II  the  inclination  of  the  interface 
between  the  two  streams  is  the  greatest  (crossing  fT/2)  reducing  to  a  local 
minimum  the  duration  of  the  interface  crossing  and  therefore  .  Neither 
the  phase-locked  streaklines  nor  the  smoke  pictures  which  are  locked  in 
phase  to  the  flap  oscillations  exhibit  any  pairing  process,  so  often  dis¬ 
cussed  in  the  literature. 

We  have  compared  streaklines  calculated  from  velocity  signals 
phase-locked  to  the  disturbence  frequency  with  velocity  signals 
phase- locked  to  its  first  subharmonid,  and  the  two  shovred  no  distinct 
difference.  The  the  first  subharmonic  of  the  oscillations'  frequency  has 
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no  influence  on  the  phase-locked  streaklines  pattern. 

Fig.  4  shows  a  sequence  of  streakline  patterns  plotted  with  an  inter¬ 
val  of  (T/2  between  then.  One  can  thus  observe  the  time-development  of  the 
flow  field,  the  changes  in  the  inclination  angles  of  the  structures  while 
moving  from  region  I  to  region  II,  and  the  stretching  of  the  structures  in 
region  III. 

The  streakline  patterns  calculated  from  the  measured  velocities  give  a 
true  descriptive  picture  of  the  phase-locked  flow  and  can  be  used  in  com¬ 
paring  flow  visualization  to  measurements  of  dynamically  significant  quan¬ 
tities  such  as  the  vorticity  field,  the  Reynolds  stresses,  and  the  power 
spectra. 


3. 2  Vorticity  field 

The  isodynes  ( iso- vorticity  contours)  presented  in  Fig.  5  were  calcu¬ 
lated  using  the  methods  described  in  section  2.3,  and  they  show  the 
phase- locked  vorticity  field  corresponding  to  the  time  (t*)  for  which 
streaklines  of  Fig.  4  were  drawn. 

Each  streakline  structure  concentrates  excess  vorticity  in  it  and, 
therefore,  may  be  loosely  referred  to  as  vortex.  There  are  different  dis¬ 
tributions  of  vorticity  in  the  three  regions  :  A  vortex  in  region  I  con¬ 
tains  two  cores,  displaced  both  laterally  and  in  the  streanwise  direction. 
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This  description  agrees  with  predictions  based  on  linear  stability  theory 
as  given  by  Michalke  0  565).  Each  of  the  vortices  seen  in  region  II  con¬ 
tains  only  one  core.  The  vortices  in  region  III  are  being  stretched  in  the 
longitudinal  directiori,  changing  their  shape  from  being  approximately  cir¬ 
cular  to  being  elliptical  with  an  ever  increasing  major  axis. 

Fig.  6  shows  a  squence  of  i so- vorticity  contours  plotted  at  TT/2  in¬ 
terval  between  them,  plotted  at  the  sane  points  as  Fig.  4  The  vorticity 
field,  as  well  as  the  streaklines'  pattern,  is  periodic,  with  a  time  period 
of  T. 

The  inclination  angle  of  a  vortex  is  related  to  the  sign  of  the  Rey¬ 
nolds  stress  as  schematiclly  sketched  below  (Browand  1980).  When  the  in¬ 
clination  angle  of  the  eddy  o(  >  n/2  the  associated  Reynolds  stress  is  posi¬ 
tive  and  when  <  H7 2  the  Reynolds  stress  changes  sign,  and  becomes  nega¬ 
tive. 

uv>0 

Xis 

The  phase-locked  vortices  in  region  I  mapped  by  the  isodynes  or  the 
streaklines,  resemble  ellipses  with  a  negative  inclination  angle.  It  is 
expected  that  the  phase-locked  Reynolds  stresses  in  this  region  would  exhi¬ 
bit  positive  values,  while  the  phase-locked  Reynolds  stresses  of  region  II 
will  go  through  zero  and  would  be  negative  in  region  III,  since  the  major 
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axes  of  the  vortices  in  these  regions  have  a  positive  inclination  angle 
with  the  flow  direction.  The  phase-locked  Reynolds  stresses  are  associated 
with  the  phase-locked  energy  production  (see  Hussain  1983)  and  it  seems 
that  while  the  phase- locked  vortices  of  region  I  would  extract  energy  from 
the  mean  flow  (negative  inclination  angle,  associated  with  positive 
phase-locked  Reynolds  stresses)  ,  they  vrould  lose  energy  to  the  mean  flow  in 
regions  II  and  III.  A  more  detailed  discussion  of  the  Reynolds  stresses  in 
the  various  flow  regions  is  given  in  section  3.3. 


The  phase-locked  vortices  in  region  III  show  no  increase  in  their  la¬ 
teral  dimension  while  their  longitudinal  dimension  increases  (stretching 
process).  This  affects  the  mean  width  of  the  mixing-layer  (or  its  manentun 
thickness)  through  the  time  integral  of  the  velocity  fluctuations  that  form 
the  basis  for  the  mean  width.  The  longer  the  vortex,  the  greater  its  dona¬ 
tion  to  the  mean  values. 


3,3  Reynolds  stress  and  power  spectra 


Since  the  Reynolds  stress  =  -  u'v'  is  linked  to  the  mean  manentun 
equation,  the  changes  in  it  must  also  reflect  respective  changes  in  the 
manentun  thickness  .  Though  the  division  of  the  flow  field  into  three 
regions  of  different  behaviour  seems  to  fit  at  first  sight,  we  have  found 
that  there  are  different  characteristics  of  the  flow  within  each  of  these 
regions,  and  we  will  discuss  the  Reynolds  stress  according  to  a  more  deta¬ 
iled  division. 


Profiles  of  the  mean  -u’ v'  product  (designated  hereafter  as  -uv  )  at 
four  character istic  stream-wise  stations  within  the  part  of  region  I  that 
exhibits  linear  growth  of  the  mixing  layer  with  increasing  X,  (  200  <  X  < 
500  urn)’,  are  plotted  in  Fig.  7,  marked  by  crosses.  The  profiles  are  posi¬ 
tive,  implying  that  energy  is  transferred  from  the  mean  flow  to  fluctua¬ 
tions  in  the  flow.  Cn  the  same  graph  we  have  plotted  profiles  of  the 
phasa-locked  -u' v'  product  (designated  hereafter  as  <-uv>)  ,  marked  by  tri¬ 
angles.  The  two  are  essentially  identical,  both  are  positive,  and  the 
difference  between  the  two  is  attributed  to  random  turbulent  fluctuations 
that  have  been  filtered  out  by  the  phase- locked  averaging.  The  mixing 
layer  is  dominated,  in  this  region,  by  the  phase-locked  fluctuations  that 
produce  phase-locked  kinetic  energy  by  the  action  of  the  phase-locked  aver¬ 
aged  Reynolds  stress  against  the  mean  strain  rate  U/  Y. 

In  Fig.  8  the  cross-spectra  graphs  at  4  lateral  location^,  in  the 
same  X-stations  are  plotted  .  The  frequency  range  of  each  graph  (and  any 
apectral  graph  presented  hereafter  )  is  between  4  and  200  Hz,  with  a  reso¬ 
lution  of  4  Hz. 

The  dominant  frequency  in  the  cross  spectra  in  this  section  of  region 
I  is  the  disturbance  frequency,  although  the  harmonic  frequency  is  mostly 
detected,  having  a  positive  peak  as  well. 

The  power  spectra  graphs  at  the  same  locations  are  presented  in  Figs. 
9  (The  power  spectra  of  the  longitudinal  component  u  )  and  10  (the  lateral 
caoponent  v).  These  graphs  are  plotted  on  semi-logarihraic  coordinates,  the 
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frequency  axis  being  the  sane  as  the  cross-spectra’s,  and  the  power  axis 
divided  into  four  orders  of  magnitude.  The  power  spectra  graphs  reflect 
the  dominance  of  the  disturbence  frequency  and  its  higher  harmonics.  The 
region  of  linear  mcnentun  thickness  growth  is  governed,  therefore,  by  the 
phase-locked  fluctuations  of  the  disturbance  frequency,  to  which  energy  is 
transferred . 


In  Fig  11  we  have  plotted  two  distributions:  the  crosses  represent 
the  local  extramum  of  each  -uv  profile,  and  the  triangles  the  local  extra- 
mun  of  the  <-uv>  profiles  vs.  the  stream-wise  coordinate  X.  These  graphs 
emphasize  the  resemblence  between  the  mean  and  phase-locked  characteristics 
in  this  sub-regiort,  and  reflect  the  difference  caused  by  turbulence. 


The  transition  region  between  region  I  and  II  500  <  X  <  560  mm)  is 

characterized  by  the  reduction  of  the  mixing-layer's  growth  rate  to  0  at 
X=560  mm.  Fig.  12(a)  represent  profiles  of  the  mean  and  phase-locked 
-u'v'  products  at  X=520  mm.  A  change  in  the  sign  of  both  profiles  can  be 
detected  at  the  upper  side  of  the  profiles  (faster  stream  side)  and  both 
are  negative  there.  Fig.  13(a)  shows  the  Reynolds  stress  spectra  at  the 
same  X-location,  and  the  disturbance  frequency  appearing  there  has  a  posi¬ 
tive  peak,  whilst  the  first  harmonic  is  negative  at  the  faster-stream  side 
and  positive  at  the  slower-strean  side.  At  X=560  mm  (Fig  12(b))  the  maxi- 
mun  positive  value  of  each  profile  at  the  faster-stream  side  is  equal  to 
the  maximum  negative  value  at  the  slower-strean  side.  At  this  location  the 
first  harmonic  exhibits  negative  peaks  all  along  the  cross-section  (Fig 
13(b)),  and  the  first  negative  peak  appears  at  the  faster-stream  side  at 


fWJ 


.  .  A  .* 


Page  IS 


the  imposed  frequency.  The  trend  of  the  disturbance  frequency  to  perform 
negative  -u'v'  product  values  continues  downstream  to  X=640  mm,  where  the 
last  peak  at  the  slower-strean  side  becomes  negative.  (Fig.  13(c)  and 
(d)).  This  trend  appears  at  both  types  of  the  -u'v*  profiles  as  a  continu¬ 
ous  transition  of  the  values  at  the  slower-strean  side  to  become  negative. 
The  energy  spectra  of  both  velocity  components  (Figs.  14  and  15)  at  the 
same  X-station  still  reflect  the  dominance  of  the  oscillations  frequency 
and  its  higher  harmonied,  over  the  flow. 

The  transition  of  the  mixing-layer  from  a  region  of  growth  to  a  region 
of  shrinking  in  size,  involves  an  interaction  between  the  mean  flow  and 
phase-locked  fluctuations  of  both  the  disturbance  frequency  and  its  first 
harmonic.  It  seems  that  the  reduction  in  the  mixing-layer's  growth-rate 
begins  at  a  point  where  the  first  harmonic's  component  of  the  -u' v'  product 
changes  sign  ,  and  becomes  negative  at  the  faster-stream  side.  The  fluctu¬ 
ations  at  this  frequency  lose  energy  to  the  mean  flow,  whilst  the 
phase-locked  fluctuations  at  the  imposed  frequency  are  still  extracting  en¬ 
ergy  from  it  (characterized  by  positive  peaks  at  the  cross-spectra  graphs). 
The  balance  between  the  two  across  the  flow  determines  the  rate  of  growth 
of  the  mixing-layer.  At  stations  where  the  components  of  the  oscillations 
frequency  are  positive  and  higher  than  the  nagative  values  of  the  first 
harmonic  components  (  500  <  X  <  560  mm),  the  phase-locked  profile  is  partly 
positive  and  partly  negative  and  the  positive  values  are  higher  than  the 
negative  absolute  values.  The  integral  of  the  profile  over  the 
cross-section  is  consequently  positive,  causing  the  mixing  layer  to  grow  at 
a  slower  rate  than  at  earlier  stages  of  region  I.  At  X=560  mm  ,  where  the 


the  imposed  frequency  turns  negative  too,  there  is  some  equilibria  between 
the  positive  phase-locked  behaviour  at  the  faster- stream  side  and  the  nega¬ 
tive  behaviour  at  the  slower- stream  side  causing  an  inhibition  in  the  mix¬ 
ing-layer’s  growth.  Oaring  the  transition  of  the  imposed  oscillations  fre¬ 
quency  toward  negative  values,  the  phase-locked  profiles  also  become  nega¬ 
tive,  and  since  the  phase-locked  phenomena  governs  the  flow  as  reflected  by 
the  resemblence  between  the  two  types  of  profiles  in  this  sub-region,  it 
causes  an  energy  transfer  from  the  phase-locked  fluctuations  to  the  mean 
flow,  followed  by  a  reduction  in  the  mixing  layer’s  width. 

Thud,  instead  of  dividing  the  mixing  layer  into  a  region  of  growth 
(region  I)  and  a  region  of  shrinking  (region  11)  we  have  divided  it  to  a 
region  of  linear  growth  (  X  <  500  nm) ,  a  transition  region  (  500  <  X  <  640 
mm),  and  a  region  of  shrinking  that  will  be  dicussed  below.  The  transition 
region  is  characterized  by  a  gradual  change  of  the  growth  rate  from  the 
constant  value  of  the  linear  region  to  the  negative  value  of  the  shrinking 
region . 

The  region  of  decreasing  manentum  thickness  (  640  <  X  <  840  mm)  is 

characterized  by  negative  profiles  of  both  -uv  and  <-uv>  (Fig.  16).  The 
absolute  value  of  each  minimun  reduces  with  increasing  X,  the  value  of  the 
mean  minimun  decreases  faster  than  the  phase-locked.  The  first  mean  posi¬ 
tive  extramun  appears  at  X=840  n*n  .  The  difference  between  the  two  types 
of  profiles  increases  with  increasing  X  and  is  attributed  to  non-locked 
fluctuations  at  frequencies  lower  than  the  imposed  frequency  that  appear  in 
the  flow  and  are  seen  at  the  cross-spectra,  u-power  spectra  and  v-power 


spectra  graphs  (Figs.  17,  18  and  19,  respectively). 

The  negative  values  of  the  -uv  profiles  in  this  region  are  in  accor¬ 
dance  with  the  shrinking  of  the  momentum  thickness,  most  of  the  energy  lost 
to  the  mean  flow  is  due  to  the  phase-locked  frequencies  (the  disturbance 
frequency  and  its  first  harmonic  )  exhibiting  negative  peaks  at  the 
cross-spectra  graphs  of  Fig.  17  that  cause  the  phase-locked  -u'v'  profile 
to  be  negative  all  over1,  but  fluctuations  at  lower  frequencies  that  are  not 
locked  in  phase  with  the  imposed  oscillations  extract  energy  from  the  mean 
flow.  Their  contribution  is  still  rather  snail. 

mm  is  another  point  of  change  in  the  mixing  layer’s  behaviour, 
and  it  marks  the  boundary  between  the  region  of  shrinking  mixing  layer  (re¬ 
gion  II)  and  the  region  of  the  regrowth  of  the  mixing  layer  linearily  with 
X  (region  III).  But  within  region  III  we  have  found  a  second  transition 
sub-region  bounded  between  8^0  <  X  <  1080  mm. 

This  transition  sub-region  is  characterized  by  positive  values  that 
appear  at  the  phase-locked  -u’v’  product  profiles  (Fig,  20).  The  mean  -uv 
profiles  remain  positive  throughout  this  region  and  downstream,  and  are  in 
accordance  with  the  growth  of  the  mixing-layer.  The  phase-locked  positive 
values  imply  that  part  of  the  large  eddies  (their  center)  is  gaining  energy 
fhom  the  mean  flow,  while  other  parts  are  losing  energy.  The  important 
role  of  the  positive  values  of  the  phase-locked  profiles  is  demonstrated  in 
Fig.  11  where  the  extrema  of  the  <-uv>  profiles  appear  to  be  positive  in 
this  regiori,  between  negative  extrema  up  and  downstrean.  This  behaviour 
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can  be  explained  in  view  of  linear  stability  considerations  that  predict 
that  whenever  the  momentum  thickness  of  the  mixing-layer  is  smaller  than  a 
given  value  (corresponding  to  the  value  of  the  neutral  case)  the  layer  is 
unstable  to  disturbances  in  the  .mposed  frequency.  Thus  it  is  possible 
that  the  vortices  of  the  oscillation's  frequency  are  gaining  energy  due  to 
the  mixing  layer's  instability,  up  to  a  point  where  the  momentum  thickness 
reaches  its  saturation  value  (  X  =  1080  mm).  The  first  harmonic  exhibits  po¬ 
sitive  peaks  at  the  Reynolds  stress  spectra  starting  at  X=860  and  downstre¬ 
am,  its  value  decreasing  with  increasing  X.  The  Reynolds  stress  spectra 
and  the  power  spectra  graphs  (Figs.  21 ,22  and  23)  exhibit  the  dominance  of 
the  disturbance  frequency,  and  the  existence  of  non- phase-locked  frequen¬ 
cies. 


The  rest  of  region  III  (  1080<  X  mm)  is  characterized  by  the  following 
features:  The  mean  -u' v'  product  profiles  are  all  positive,  and  in  accor¬ 

dance  with  the  mix ing- layer '  s  growth  (Fig.  21*).  The  <-uv>  profiles  are 
negative  and  imply  that  the  phase-locked  vortices  are  losing  energy  to  the 
turbulent  mean  flow.  They  do  dominate  the  Reynolds  stress  spectra  (Fig. 
25) ,  but  other  non  phase-locked  frequencies  that  exhibit  lower  peaks  appear 
and  are  responsible  for  the  mixing-layer's  growth.  It  may  be  noted  that 
the  two  distinct  frequencies  that  appear  in  this  region  are  12  and  3 2  Hz. 
Nowhere  along  this  region  does  the  first  subharmonic  of  the  imposed  fre¬ 
quency  appear.  Since  this  frequency  is  assuned  to  govern  phase-locked  eddy 
amalganationS,  it  is  obvious  that  such  a  mechanism  would  not  exist  in  this 
mixing-layer.  The  regrowth  of  the  mixing-layer  in  region  III  can  be  attri¬ 
buted  to  either  random  amalgamations  of  vortices,  or  some  non-linear  inter- 
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action  of  fluctuations  at  the  lower  frequencies,  that  are  not  locked  in 
phase  with  the  disturbance.  This  conclusion  is  in  agreement  with  the  con¬ 
clusion  of  Ho  and  Huang  (1980),  that  the  subharmonic  is  the  agent  rather 


than  the  product  of  the  vortex  amalgamations. 
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F igure  Captions 


Particle  Streaklines. 

Comparison  between  Calculated  Streaklines  and  Smoke  Picture 
Three  Regions  of  the  Mixing  Layer. 

Sequence  of  Streaklines. 

Iso-Vorticity  Contours. 

Iso-Vorticity  Contours. 

Mean  uv  product. 

Q*oss  Spectra  Distributions. 

Power  Spectra  Distributions  of  U. 

Power  Spectra  Distributions  of  V. 

Distributions  of  -uv  and  <-uv>. 

Mean  and  Jhase  Locked  -uv  Products. 

Reynolds  Stress  Spectra. 

Energy  Spectra  of  U. 

Ehergy  Spectra  of  V. 

Momentum  Thickness. 

Cross  Spectra. 

U  Power  Spectra. 

V  Power  Spectra. 

Transition  Sub  Region. 

Reynolds  Stress  Spectra. 


U  Power  Spectra 
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2 3.  V  Rower  Spectra. 

24.  Mixing  Layer  at  Region  III. 

25.  Reynolds  Stress  Spectra. 
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